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ABSTRACT 


This thesis deals with analytical predictions and experimental determination of the 
modal frequencies and shapes associated with the arm portion of the Naval Postgraduate 
School's Flexible Spacecraft Simulator (FSS). A description of how piezoceramic sensors 
and actuators are incorporated in finite element modeling is presented. A MATLAB™ 
code conducting the finite element modeling of the arm is used to verify the modal 
frequencies generated via Structural Dynamic Research Corporation's I-DEAS™ 
software. Modal testing is conducted with the I-DEAS Test package to determme the 
first four physical modal frequencies and shapes so that we may compare them with the 
analytical results. The results of the testing indicated that finite element analysis predicted 
modes one and three within an average of 23.1%, and modes two and four within an 


average of 4.4%. 
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I, INTRODUCTION 


The often times prohibitive cost for putting payloads into orbit has lead the 
scientific and engineering community to design lightweight yet dynamically stable 
structures. Accurately predicting the natural frequencies of a structure, corresponding 
mode shapes and modal damping is integral to the design and dynamic control of such 
structures. Lightweight, flexible structures inherently possess low resonant frequencies 
which can easily be stimulated by a variety of spaceborne excitation sources such as 
attitude control thrusters and reaction wheels. 

Previous work done on flexible structural control utilizing Positive Position 
Feedback (PPF) and Strain Rate Feedback (SRF) by Newman and Feuerstem, [Refs. 1 
and 2] mandated an accurate knowledge of modal characteristics for successful control 
law design. This previous work in modal control provides the stimulation for this 
research. 

The most common method of modal frequency and shape prediction is some form 
of the Finite Element Analysis (FEA) approach. This method exploits digital computing 
power and provides a formulation for the analysis of dynamic characteristics of flexible 
structures. By regarding a complex structure as a composite of finite elements, the 
method requires that displacements be compatible and internal forces in equilibrium at 
points known as nodes, which are shared by several elements. The result being the total 
structure should act as if one continuous entity. 

The Naval Postgraduate School Flexible Spacecraft Simulator's arm provides a 
testbed for modal analysis. It provides a reasonably accurate representation of a very 
flexible arm such as that which might be utilized for a deployed antenna, possibly imposing 
very accurate pointing constraints. This arm employs nearly collocated piezoceramic 
actuators and sensors for the damping and control of vibration. The FEA method is well 


documented, [Refs. 3 and 4]. What is initially explored is how the physical properties and 


forces implied to a structure via the incorporation of these piezoceramic sensors and 
actuators modifies the FEA. 

Many commercial software FEA packages are available to assist with the 
predictions described above. For the purpose of this thesis, the modal prediction package 
provided by Structural Dynamic Research Corporation's I-DEAS software is used. A 
comparison of this output to a simple MATLAB code which utilizes standard FEA 
techniques [Ref. 4 ] is then undertaken. 

Finally, modal analysis of the arm is conducted experimentally so that we may 
determine the accuracy of analytical predictions. This modal testing utilizes the testing 
package incorporated within the I-DEAS software. The arm is excited at several nodes 
with an impact hammer with a built-in force transducer. Accelerations of the arm are 
measured via an accelerometer whose signal has been conditioned and amplified. These 
analog time domain signals are discretized and then changed into the frequency domain via 
the discrete Fourier transform. The analyzed signals can then be manipulated m a variety 
of ways to produce such information as natural frequencies, mode shapes and damping 


ratios. 





I. FINITE ELEMENT ANALYSIS WITH PIEZOCERAMIC 
ACTUATORS AND SENSORS 


A. ANALYTICAL BACKGROUND AND ASSUMPTIONS 
For the purpose of discussing the procedure with which FEA methods are 


modified to take into effect the physical and active characteristics of an element 
incorporating piezoceramics, consider the co-located sensor/actuator pair on the n-th 
diemest of an isotropic flexible beam as depicted in Figure 1 below. All physical and 
electrical properties and dimensions of the actuator and sensor are assumed identical. 

The neutral axis is aligned with the centerline of the beam. The beam is free to translate 


in the +/- number 3 axis direction and rotate about the number 2 axis only. 


node n-1 —_» »< 


n-th element 
co-located actuator and sensor 





Figure 1. Piezoceramic Actuator and Sensor, n-th Element 


B. GOVERNING EQUATIONS 
1. Passive Contributions 

We first consider a structural element as shown in Figure 2. The equations of 
motion for the element are of the form: 

Mq+Kq=0 (1) 

where V/ and K are the elemental mass and stiffness matrices, respectively. q represents 


an element displacement vector of the form: 


Wn-l 


Ort (2) 


where W,-) and W,, are the transverse deflections of the left and right end of element n, 
respectively, and @,-; and 9, are the rotations of the left and right end of element n, 


respectively. These displacements are shown in Figure 2 below. 






piezoceramic 


on ats 


theta piezoceramic \ theta 
n-| Nl 









Ww Ww 
n-l n 


k ———| 


Figure 2. n-th Element Displacements 





By deriving the equations of motion for the element using Lagrange's equations, which 
reduces, in turn, to determining the kinetic energy, the potential energy, and the virtual 


work expressions in terms of the nodal cordinates, consistent element mass and stifthess 


matrices can be generated. 
The consistent element mass and stiffness matrices are given by the following: 








156 22h 54 -13h 
22h 4h? 13h -3h? 
M=2 3 
20) 54 13h 156 —22h " 


13h —3h? -22h 4h? 


12 6h -12 6h 

6h 4h? -6h 2h? 
K=< 4 
h’) -12 -6h 12 -6h “) 


6h 2h? —6h 4h? 


where m is the linear mass density of the element, 1 is element length, / is Young's 
modulus for the element and / is the element moment of inertia about the bending axis. 


[Ref. 4 | 
In order to account for piezo contribution to the mass matrix, ™ must be equal to 


the linear mass density of the beam plus the linear mass density of the piezoceramic 


actuator and sensor, or, 


mM = Mbeam + M piezo (S) 
To account for piezo contribution to the stiffness matrix, the // term must be of the form: 
EI = Eoecaml beam + E piezol piezo (6) 
Since the neutral axis is aligned with the beam centerline, [4 cam is simply: 
wt? 
Ibeam =z (about bending axis) (7) 


where W is the width of the element and ¢, is the thickness of the structure. 
By employing the parallel axis theorem to determine [piezo, we ative at the 
following: 
wt th t 
Ipiezo = Gy + Wty +z)" (8) 
where f,, is the thickness of the piezoceramic. 
After simplification, 


. t 
lies =i, & +4 : (9) 


Eeam and E piezo are simply the Young's modulus values for the beam and piezo, 


respectively. 


Once the particular elements incorporating the piezos are modified as above, all 
other aspects of constructing the global mass and stiffness matrix for an entire structure 
incorporating multiple elements are straightforward [Ref 4]. The global version of 
Equation (1) will now present an eigenvalue problem whose solution will reveal modal 
frequencies and shapes. 

Ze Active Contributions 

a. Sensor 
Piezoceramic actuator and sensor theory is well documented [Ref 5]. 


Voltage output from the piezoceramic sensor is given by the following equation: 


Vos t (Epd31/D)(E1 + €2) (10) 
where 


tp is the thickness of the piezoceramic sensor 

E,, is Young's modulus for the sensor in N/m? 

di) is the lateral strain coefficient in m/V or Coul/N 
Dis Abs Permittivity in Farad/m or N/V? 

€ 1 and €2 are longitudinal and lateral strain, respectively 


Since longitudinal strain €, and the lateral strain €2 are related by Poison's ratio, V, given 
by: 
€1V =-—€2 (11) 
Equation (10) may be written as: 
y= tp)(E,d31/D)C0 _ V)E} (12) 
Figure 3 illustrates the sensor mode of the piezoceramic. The illustration on the left of 
Figure 3 depicts the case when the sensor is subject to lateral tension and vertical 
compression. The right illustration shows the sensor in vertical tension and lateral 


compression. Produced voltage polarities under these conditions for the mdicated poling 


direction are as shown. 








poling 
gn voltage on 
F polarity F 


=| before external force 


|_| after external force 





Figure 3. Piezoceramic Sensor 


In generalized coordinates, sensor voltage produced at the n-th element is 


represented by the following, [Ref. 5]: 
Vi=V 0101 |g (13) 


where )’, is given by Equation (12) and the displacement vector q, is given by equation 
(2). 

b. Actuator 

The equation with which the general behavior of the element, exposed to 
external excitation forces and also the forces imposed by the actuator, is represented by 


the following: 
Mq + Kq = F external + Pactuator (14) 


where Factuator is given by the following expression, [Ref. 5]: 
F actuator — —B Va (15) 


V, is the voltage applied to the actuator and B is a 4 x 1 vector possessing moment 


correction terms whose units are in Newton- meters/volt: 


0 
by 
n= 16 
0 (16) 
b, 
with 
b> = —d3Epw(2 + 2 and b4= b> (17) 


where all variables are as in Equation (10) and w is the actuator width. Figure 4 
illustrates the actuator mode of the piezoceramic. The left side illustration demonstrates 


vertical contraction and the corrective moment producing lateral expansion when applied 
voltage is of opposite polarity from the poling voltage. The right side demonstrates 
vertical expansion and the corrective moment producing lateral contraction when a 


voltage of the same polarity as the poling voltage is applied to the piezo. 








applied voltage applied voltage 
opposite polarity same polarity 
of poling voltage of poling voltage 


V- Vt 
: 


poling 
voltage 
polarity 


gn 


Ss before applied voltage 


[| after applied voltage 





Figure 4. Piezoceramic Actuator 


C. Coupling Sensor Output to Actuator Input 

The last procedure is to determine the control scheme which will couple 
the measured sensor voltage output to that applied to the actuator in order to produce a 
restoring moment to the element and ultimately the entire structure itself. We will 


examine a simple implementation of velocity feedback, where 

V,=kV" (18) 
k is a voltage gain factor. Combining all variables in Equations (2), (13) and (15) will 
ultimately lead to the following overall expression for Factuator? 


v Wr] 
b Z 6,1 

F actuator = —k 0 V.| 0101 | Ws (19) 
b4 0, 


With the variables 52 , b4 and V; as before. Equation (19) can now be written in the 


form: 
P actuator = —Cq (20) 
where C is a 4 x 4 matrix given by: 
| v 
by 
C=-k : 21 
5 (ML O101 | (21) 
b4 


With no external excitation forces present, Equation (14) may now be rewritten mm 


the general damped form: 
Mg+Cq+Kq=0 (22) 


In the construction of a global system of i number ofelements, the global C 


matrix will be null in all segments except those which possess actuators. Therefore, the 
passive damping of the structure is neglected. Equation (22) then presents a 
conventional eigenvalue problem which will reveal modal frequencies, modal vectors and 


damping ratios. [Ref. 6] 
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I. FINITE ELEMENT ANALYSIS COMPUTER PREDICTIONS 


A. I-DEAS™ FEA MODELING 

I-DEAS is a menu driven, Windows™ type, software package developed by 
Structural Data Research Corporation. Its finite element analysis and testing packages are 
just two of several structural design tools available to designers and engineers. 

The first step in conducting FEA with I-DEAS is to model the structure under 
consideration. Figure 5 depicts the mesh geometry for the Flexible Spacecraft Simulator's 
arm used for the analysis. This flexible arm rests on a granite table which simulates a 
friction free environment and is supported at the elbow and tip by two airpads (nodes 10 
and 19). Node 1 is rigidly fixed to the table. Mesh geometry is such that the arm is free 
to translate in the +/- number 1 and 3 axis and rotate about the number 2 axis only. 

The arm is constructed of 7075 T-6 aluminum, with a cross section of 1/16" x 1" 


(ty x W), and has several mass intensifiers, which are treated as pomt masses for the 
purpose of the analysis, spaced throughout the beam. Elements 1, 10, 11, and 18 have 


two wafers per side of Navy Type II PZT bonded to the beam surface. Table 1 below 


cites the material and electrical properties of this piezoceramic. 


Description Units Value 
Lateral strain coeff m/V or Coul/N 1.8 e-10 


Young's modulus Pascal 6.3 e10 
Poison's ratio n/a 0.35 

Abs Permittivity Farad/m or N/V? 1.58 
piezo thickness mM 1.905 e-4 





Mass density kg/m? 7.7e3 


Table 1. Navy Type II PZT Properties 


i] 


ARM STRUCTURE ON GRANITE TABLE 


n = nodes 
(Airpads at nodes 10 and 19) e = elements 
n-1 n2 n3 n4 n5 n-6 n-/7 n-8~) 1-9 n-10 
e-] e-2 7 e-3 -—ye-4 -e-5 —e-8 ee-7 mo e-8 —e-9 
| e-10 
n-11 
BEAM MATL: 7075 T-6 e-11 
| - CROSS SECTION: 1/16 " x 1", (tb x w) 12 
- DENSITY: 0.003144 slugs/in3 wt 
-E=1.029E7 psi 
e-12 


MASS INTENSIFIERS: [| = 0.0317 sugs 


AIRPAD: | = 0.0252 slugs 
AIRPAD: = = 0.0567 slugs 


NODAL COORD (1,3 AXIS) in.: 


n-1 (0,0) n-11 (26.5,-2.5) 
n-2 (4.75,0) n-12 = (26.5,-5.0) 
n-3  (7.25,0) n-13 = (26.5, -9.87) 
n-4 (9.75,0) n-14 = (26.5,-12.37) 
n-5 (12.25,0) n-15  — (26.5,-14.87) 
n-6  (14.75,0) n-16 = (26.5,-17.37) 
n-7  (17.25,0) n-17  — (26.5,-19.87) 
n-8 (19.75,0) n-18 = (26.5,-21.5) 
n-9  (24.0,0) n-19  (26.5,-26.5) 


n-10 (26.5,0) 
mm ACCELEROMETER =0.005 slugs 


NAVY TYPE II PZT ON 
ELEMENTS | e-l, e-10, e-11, e-18 





Figure 5. Finite Element Analysis Mesh for 


Flexible Spacecraft Simulator Arm 
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The mesh geometry (all nodal coordinates), beam element geometry, element 
material properties, point mass locations and values, and structural restraints are all mput 
into the I-DEAS modeler. Once this model is stored, modal analysis can be conducted. 
I-DEAS constructs element mass and stiffness matrices, assembles the global mass and 
stiffress matrices for the structure and then calculates the modal frequencies and 
associated mode shapes. Modal analysis was conducted on models with and without the 
PZT wafers on elements 1, 10, 11, and 18. Appendix A contains the -DEAS printouts 
for the model, and the deformed geometry for the first four modes, with and without 
piezoceramic contribution on elements 1, 10, 11, and 18. Table 2 below lists the 


frequencies associated with the first four modes generated by I-DEAS FEA for the arm. 


Frequency (Hz) Frequency (Hz) 
(without piezos) (with piezos) 


0.136 0.159 
0.372 0.469 
1.651 1.87 
2.373 3.06 





Table 2. I- DEAS FEA Modal Frequencies 


B. MATLAB™ FEA MODELING 

The MATLAB code used to conduct FEA on the FSS arm appears m Appendix 
B. No control laws were implemented relating sensor output to actuator mput. The 
code, therefore, performs the eigenvalue problem shown in Equation (2). This code was 
also run with and without the PZT properties included in order to compare it with the 
J-DEAS output. Appendix C contains printouts of the modal frequencies, with and 
without piezos, in descending order. Table 3 below lists the modal frequency predictions, 
for the first four modes, generated from MATLAB. The I-DEAS predictions are also 


re-listed for comparisons sake. 
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MATLAB I-DEAS MATLAB I-DEAS 
(with piezos) (with piezos) (without (without 
piezos) 

(Hz) (Hz) 





Table 3. MATLAB and I-DEAS FEA Modal Frequencies 


As seen from the table) MATLAB and I-DEAS modal frequency predictions for 
both cases correlate well. For the case with piezo contribution to the arm considered, the 
I-DEAS predictions diverge slightly from MATLAB. The MATLAB and I-DEAS 
predictions for the arm without piezo contribution considered are essentially the same. 
The overall effect of the piezos on the arm stiffens the structure incrementally, as 


evidenced by the slightly higher frequencies generated by both MATLAB and I-DEAS. 
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IV. MODAL TESTING 


A. INTRODUCTION 

In Chapters II and HI, the differential equations, or models, of the system are 
assumed to be known, and the theory developed consists of calculating and characterizing 
the response of the system to known inputs. This is called the "forward problem". In this 
chapter the interest lies in measuring the response of a structure and m some way 
determining the equations of motion from the test data. The problem of determming a 
system of equations from information about inputs and responses belongs to a class of 
problems called "inverse problems", or system identification. The measurement of the 
vibrational behavior of the structure is used to verify the mathematical model of the test 
structure outlined in Chapters I and IIE. [Ref 3] 

Modal testing is the interpretation of test data collected from a vibrating structure. 
The purpose of a modal test is to construct a mathematical model of the vibrational 
properties and dynamic behavior of the tested structure. 

Modal testing focuses on frequency, shape determination and dampmg. On a plot 
of frequency response as a function of frequency, a resonance frequency (modal 
frequency) is marked by a peak of the magnitude. A second phenomena of resonance 1s 
that the phase of the response shifts by 180 degrees as the frequency sweeps through 
resonance, with the value of the phase at resonance bemg 90 degrees. 

B. MEASUREMENT HARDWARE 

Vibration measurement generally requires several hardware components. These 
consist of a source of excitation for providing a known and controlled mput to the 
structure, a transducer or accelerometer which converts mechanical motion of the 
structure to an electrical signal, and a signal conditioning amplifier to match the signal 
level of the accelerometer to that of the data acquisition system. This data acquisition 
system provides the interface between excitation forces and accelerometer signals to the 


testing software imbedded within I-DEAS. 
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Appendix D delineates all hardware associated with modal testing of the FSS arm. 
An impact hammer with a built-in force transducer was used to excite the structure. 
Although the impact hammer is simple and adds no mass to the arm, it sometimes may be 
incapable of transforming sufficient energy to the structure to develop adequate response 
signals in the frequency range of interest. | Consistent impact techniques can also be 
difficult to develop. The key frequency range of interest for this test was from 0 to 4 Hz. 
This bandwidth would hopefully contain resonance information for the first four flexible 
modes. 
C. DIGITAL SIGNAL PROCESSING 

Modal Testing analysis is performed in the frequency domain. The data acquisition 
system converts the analog time domain signals into the digital frequency domain via the 
Discrete Fourier transform. This data is then passed on to I-DEAS in order to perform 
the required computations. Briefly, by comparing the transfer functions of the response of 
the structure to the transfer function of the excitation force, one can determine the transfer 
function of the structure itself and ultimately its associated resonances. Mode shapes may 


also be extracted from this data. The transfer function relationship is given by: 


G(j@) = Sa. V@V/Sp(V) (23) 
where, 
Sy(/®)= response of the structure in the frequency domain, @ 
S(i@)= excitation force in the frequency domain, @ 


G(j@)= the transfer function of the structure 


Reference 7 provides an excellent source for information on the theoretical and practical 


problems surrounding the measurement, processing, and analysis of modal test data. 
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D. TEST PROCEDURES 

A modal test mesh must be created within the I-DEAS testing environment. This 
mesh is depicted in Figure 6. Unlike the mesh created for the finite element analysis, no 
material properties need be entered. The nodal coordinates in this mesh defines the 
degrees of freedom for the structure. The accelerometer was located at nodes 13, 10 and 
8 for the bulk of the testing due to physical constraints (proximity to the piezos and 
associated wiring). Node 16 was used as the primary impact location for defining the first 
four modal frequencies. Impact at this location proved to consistently and reliably excite 
the first four modal frequencies. Impact at all other nodes was performed in order to 


determine mode shapes. [Ref. 8] 
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(26.5,-17.37) 
(26.5,-19.87) 
(26.5,-26.5) 


Figure 6. Modal Test Mesh for 


Flexible Spacecraft Simulator Arm 
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E. TEST RESULTS 

Figures 7, 8 and 9 are I-DEAS printouts of the arm's frequency response function 
measured at nodes 13, 10 and 8 respectively. Impact for all three of these measurements 
took place at node 16. The resonant frequencies are clearly represented by distinct peaks. 
The 180 degree phase shift (or phase reversal) crossing resonance is also reasonably well 
defined. Table 4 below lists the averages of the modal frequencies observed during these 
tests at nodes 13, 10 and 8. The table also re-lists the MATLAB FEA predictions (with 


piezoceramic contribution) and a percent deviation. 


Modal Test MATLAB FEA % 
Frequency Averages Prediction Deviation 
(impact at node 16, (with piezos) 
response at 
nodes 13, 10 and 8) 


(Hz) 





Table 4. Modal Test Averages and % Deviation from FEA 


Figure 10 depicts an analytical curvefit that was performed by I-DEAS on the test 
data obtained at node 13. This analytical curve is distinguished from the raw data by its 
smooth character. The peaks identified by this curvefit are essentially identical to those 
shown on Figure 7. The frequency data extracted from this curvefit was input to the mode 
shape subroutine within the I-DEAS Testing package. 

Figure 11 compares the third mode shape generated by I-DEAS FEA to that 
actually observed during modal testing. The bottom illustration of Figure 11 is the third 
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Figure 8. Modal Test Response Function (impact at node 16, accelerometer at node 10) 
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Figure 10. Modal Test Response Function Analytical Curvefit 
(impact at node 16, accelerometer at node 13) 
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ChLexacrn 
LOAD SET: 3} WOGE: J FREQ: 1.64914037 
DISPLACEMENT - RORMAL WLM: 8.00 WAR: $9.42 


FEA model 


MODER: 1 FREQ: 1.2497907 
OLSPLACEKENT - NORMAL BHINt 9.990 MAX: 1360.30 


Modal Test 





Figure 11. Mode 3 Shape (FEA and Modal Test) 
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mode shape generated during the testing, the top represents the prediction of mode 3 
behavior via FEA. 
A. structure's inherent damping characteristics can only be determined 


experimentally. An n degree of freedom system has n damping ratios, G;, 


(dimensionless) one associated with each modal frequency, @;, (radians per second). The 


free vibration characteristics of each mode decays over time, t, according to the term: 
eit (23) 


[Ref. 9]. Table 5 lists the damping ratios for the first through fourth modes obtained 
from the I-DEAS testing package. The exponential time decay coefficient is also listed. 


@; = 2n(freq) 


(radians/s) 





Table 5. Modal Test Damping Ratios 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The major reason for conducting model testing is for comparison and verification 
of analytical models. What is ultimately desired is a mathematical model of the structure 
under consideration for the purpose of predicting how the structure will behave under a 
variety of different loadings, to provide the plant in a control system design, such as 
Positive Position Feedback (PPF) or Strain Rate Feedback (SRF), and to aid m the design 
process in general. 

Modal testing of the FSS arm provided valuable insight into this field. Testing 
results yielded resonant frequencies that deviated from prediction by as little as 3.5% for 
mode four, and as much as 23.4% for mode one. Overall, modes one and three deviated 
from prediction on the average of 23.1%, and modes two and four deviated on the 
average of 4.4%. The difference between the analytical prediction and experimental 
determination of frequency for the first mode is significant. Further work is required to 
improve the model and/or improve the modal testing technique. 


The damping ratios obtained during modal testing are that of an underdamped 
system, (C; <1) , [Refs 3 and 4]. The expected dominance of mode one is supported 


by the lowest exponential time decay coefficient. 

As mentioned earlier, the most common application of modal testing is the 
comparison of measured vibration modes with those predicted by a finite element or other 
theoretical models. Resonance frequencies must be accurately identified, also a complete 
and accurate description of mode shapes is desired, especially for non-collocated sensor 
and actuator systems. Frequently, the analytical model will not predict the measured 
frequencies. As a result the analytical model may then be updated, iteratively, until it 
produces the measured natural frequencies. The modified model is then considered an 
improvement over the previous model. [Ref. 3] 

One possible model modification would be to consider the inertial contributions of 
the mass intensifiers and airpads. These were treated as point masses with zero moment 


of inertia. This would be a good starting pomt in model evolvement. 


pis 


Normally, the FSS arm is attached to the center body of the FSS which houses the 
reaction wheel and thrusters [Ref. 2]. With at least two additional "zero-Hz operation" 
accelerometers, such as that used for this experiment, placed on the arm ( total of three 
on the arm), and with an additional and suitable accelerometer mounted to the center body 
for excitation source measurement, more reliable frequency measurements may be 
expected. These results may then be used to "tweak" any implementation of PPF or SRF 


that may be used for active control, in more of a real time fashion. 
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DISPLACEMENT 





FREQ: 0. 


i-3 392-7 °5,6 
a ae | 


27 


flexarm 
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flexarm 
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$ FINITE ELEMENT 


clear 
m=1.6375e-5; 
mp=2.162e-5; 


smp=0; 


h=2.95; 
hl=4.75; 
h8=4.25; 
h12=4.8725; 
h17=1.625; 
h18=5.0; 


T=2.0345e-5; 


E=1.0298e7; 


Ip=4.56e-5; 
Ep=9.137e6; 


oP oP? 


fe) 


ae ae oo oo 


oP oP 


APPENDIX B 


ANALYSIS OF FSS ARM 


Linear Mass Density of Beam in slugs/12/in 
Linear Mass Density of Piezo in slugs/12/in 


no piezo option 


beam: 


an? 4 
1B7 ine Z 


piezo: 


ind 
Moy gal shay 


% Mass Intensifiers in slugs/12: 


Lma=0.00264; 


Lmb=0.0021; 


Lmc=0.00472; 
Accel=0.000416; 


n=18; 


# of elements 


% elemental mass and stiffness 


mMliem*h=[156. 22*h;22*h 4*h*2)/420; 
m22=m* he 1156. =22*h- 227. 24h 27 4Z0; 
mi2=m*h* (54 -13*h;13*h -3*h%2)]/420; 


m21i=m12'; 


m22e1=(m+mp) *hi* [156 -22*h1;-22*hl 4*h1*2)/420; 
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mlle8=m*h8* [156 22*h8;22*h8 4*h8%2]/420; 
m22e8=m*h8* [156 -22*h8;-22*h8 4*h8%2]/420; 
ml 2e8=m*h8* [54 -13*h8;13*h8 -3*h8%2]/420; 
m2le8=ml12e8'; 


mllell=(mtmp) *h* [156 22*h;22*h 4*h*2]/420; 
m22e11=(mtmp) *h* {156 -22*h;-22*h 4*h*2]/420; 
mi2e11=(m+mp) *h* [54 -13*h;13*h —3*h*2]/420; 
m2lell=mi2el11'; 


mlle12=m*h12* [156 22*h12;22*h12 4*h12%2]/420; 
m22e12=m*h12* [156 -22*h12;-22*h12 4*h12%2]/420; 
m12e12=m*h1i2* [54 -13*h12;13*h12 -3*h12% 21/420; 
m21e12=m12e12'; 


mllel7=m*h17* [156 22*h17;22*h17 4*h17%2)/420; 
m22e17=m*h17* [156 -22*h17;-22*h17 4*h17%2)/420; 
ml2e17=m*h17* [54 -13*h17;13*h17 -3*h17%2]/420; 
m2lel7=m12e17'; 


ml1e18=(m+mp) *h18*[156 22*h18;22*h18 4*h18%2]/420; 
m22e18=(m+mp) *h18* [156 -22*h18;-22*h18 4*h18%2]/420; 
m12e18=(m+mp) *h18* [54 -13*h18;13*h18 -3*h18%2]/420; 
m21le18=mi2e18'; 


mov=zeros (4); 
mov (1:2,1:2)=m22; 
mov (2,2)=m22(2,2)+mliell (2,2); 
mov (2,3)=m12el11 (2,1); 

nal 2)=m12e11 (2,1) 

mov (2, 4)=m12e11 (2,2); 

mov (4,2)=m12el11 (2,2) 

Hows. 4,3:4)=m22ell1; 


k11=E*I*[12 6*h;6*h 4*h*2]/h*3 
k22=E*I* [12 -6*h;-6*h 4*h*2]/h*3; 
kK12=E*I*[-12 6*h;-6*h 2*h*2]/h%3 
k21=k12'; 
kp1=Ep*Ip*[12 6*hl -12 6*h1; 

6*hl 4*h1*2 -6*h1 2*h1*%2; 

-12 -6*hl 12 -6*hi; 

6*hl 2*h1*2 -6*h1l 4*h1*2]/h1*°3; 


kp2=Ep*Ip*[12 6*h -12 6*h; 
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rRA.M 


6*h 42*h*2 =—6*h. 242; 
212 =O" 12 Pon; 
6*h 2*h*2 -6*h 4*h%2)/h%3; 


$ no piezo option 
$kpl=zeros (4); 
$kp2=zeros (4); 


k22e1=(E*1* (12. —¢6*hil;—6*hL 4*h1°2)/hi"3)+kpl (374,374); 


klle8=E*I* [12 6*h8;6*h8 4*h8%2]/n8%3; 
k22e8=E*I*[12 -6*h8;-6*h8 4*h8%2]/h8%3; 
k12e8=E*1I*[(-12 6*h8;-6*h8 2*nh8*2] /h8%*3; 
k21e8=k12e8'; 


kiiell=(B* 1" (12: 6*hs6"h 4*h°2 | 7h" 3) +kp2(122, 122); 
k22e11=(E*1I* (12 -6*h;-6*h 4*h%*%2]/h%3) +kp2 (3:4,3:4) ; 
k12e11=(E*I*[-12 6*h;-6*h 2*h%2)/h%3)+kp2(1:2,3: 4); 
k2lell=kl2ell1'; 


kl1le12=E*I* [12 6*h12;6*h12 4*h12%2)/h12%3; 
k22e12=E*I* (12 -6*h12;-6*h12 4*h12%2)/h12%3; 
k12e12=E*I* ([-12 6*h12;-6*h12 2*h12%2]/h12%3; 
k21e12=k12e12'; 


kllel7=E*I*[12 6*h17;6*h17 4*h17%2)/h17%3; 
k22e17=E*I* (12 -6*h17;-6*h17 4*h17%2)/h17%3; 
k12e17=E*I* [-12 6*h17;-6*h17 2*h17%2)/h17%3; 
k21e17=ki2el7'; 


k1lle18=(E*I*(12 6*h18;6*h18 4*h18%2]/h18%3)+kp1(1:2,1:2); 
k22e18=(E*I*(12 -6*h18;-6*h18 4*h18%2] /h18%3) +kpl1 (3:4,3:4); 
k12e18=(E*I*[-12 6*h18;-6*h18 2*h18%2]/h18%3)+kp1(1:2,3:4); 
k21e18=k12e18'; 


kov=zeros (4) ; 
kov(1:2,1:2)=k22; 
kov (2,2) =k22(2,2)+kllell (2,2); 


kov (2,3) =k12e11(2,1); 
kov (3,2)=k12e11 (2,1); 
kov (2, 4) =k12e11 (2,2); 
kov (4,2) =k12e11 (2,2); 


kov (3:4,3:4)=k22el11; 
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LOS. Lt 


& global mass matrix 


gm(1:2,1:2)=m22el+mll1; 
gom(1:2,3:4)=m12; 


for i=l:n-2 
(2*14+1:2*1+2,2*1i-1:2*1)=m21; 


(2¥i4122* 142,27 14127142) =mZ2tmil 1; 
(2*14+1:2*14+2,2*1+3:2*1+4) =m12Z; 


EEE: 


2UgL 


145138 
Oyo 
P14 45% 
-16 1.5% 


'20,19% 
rat rm Uc 
P20 Ze 
P22 pe 


724,21; 
S227 25% 
aa ae Ae ie 


PS25515 
ee TS ie 
sie (AP ee 
*34,55. 


7207302 
‘345: 35" 
$360,507 


20) 
20) 
°22) 
ZZ) 


20)=mov; 


14) =m22+mlies8; 
14)=m21e8; 
16)=m12e8; 
16)=m22e8+m11; 
=m22ell+milell; 
=m21lel1l1; 


=m22eli+milel2; 


22)=m21e12; 
24)=m12e12; 
24)=m22e12+m11; 


32)=m22+ml1llel7; 
32)=m21el17; 
34)=m12el17; 
34)=m22e17+mllels8; 


34)=m21e18; 
36)=m12e18; 
36)=m22e18; 
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FEA.M 


Lmm(11,11)=Lma; 
Lmm(13,13)=Lma; 

Lmm (23,23) =Lma; 

Lmm (25,25) =Lma; 

Lmm (29,29) =Lma; 
Lmm(31,31)=Lma; 
Lmm(27,27)=Accel; 
Lmm(17,17)=Lmb+0.0154; 


Lmm (35,35) =Lmc; 
gmm=gm+Lmm; 


% global stiffness matrix 


gk(1:2,1:2)=k22el1+k11; 
gk (1:2,3:4)=k12; 


for 1=l:n-2 
gk (2* itis 2*142,2* 14122" 1 SHK21; 


Get 2* leat 2 oe 2) SK ZZ EK 
Gk( 2° 4asLi2*i42)2* 1435227174) =KizZ; 








207.45] 


©14213% 
716,13: 
-14-15% 
164452 


220,19: 
S22, 49% 
2203-215 
i2zyp2. 


‘24,721 
S22 27235 
524523. 


See agi i lg 
34y oO) 
32y50% 
Soa y 332 


SO; 30° 


°-20)=kov; 


) 
14)=k21e8; 
16)=ki2e8; 
16)=k22e8+k1l; 
20)=k22e11+kilell; 
20)=k21lell1; 
22)=ki2el11; 


°22)=k22e1l1l+kllel2; 


>22)=k21lel1l2; 


24)=kl2el2; 
24)=k22e12+k11; 


=k22+kllel7; 


=ki2el7; 
=k?22el7+kllel8; 


32) 
:32)=k21el7; 

34) 

34) 


34)=k21e18; 
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PRA. M 


gk (33: 34,35: 36)=kizels; 
gk(35:36,35'36)=k2Ze13;7 


% mode shapes and natural freqs 


B=inv (gmm) *gk; 

[x,D]=e1ig(B); 

w=sqrt (D); 

f=w/2/pi; 

[k1,k2]=sort(diag(D)); 

x=x(:3,K2); 

% modal masses 

y=x '*gmm* xX; 

& mode shapes normalized wrt mass 
for L=l:27n 


z(:,i)=x(:,i)/sqrt(y(i,1)); 


end 
Z; 
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» fea 
» diag (f£) 


ans = APPENDIX C 
1.Q0e+003 * 


.§9118474609514 Modal Frequencies (without pi ing 
.99568135727380 qt ( t piezos), descending order 
.59480701359855 
.99704503076688 
.93678510562751 
.28716854586410 
.12001248476245 
.00095892851113 
.64377826366426 
.57337814606062 
.36667629704302 
.24571931505966 
.13139307876285 
.08173078771907 
.83137661757685 
.80808968513644 
.67562468034393 
.65960922618500 
.33282677547923 
.22223482048505 
.20688050862119 
.18556683346412 
.11445097218405 
.09656670691812 
.06764056549358 
.04693032740926 
.03658424100849 
.02811923523049 : 
.02120575875813 
.01459733742682 
.00694405426306 
.00599862204817 
.00237295384752 
.00164563740851 
.00037335596675 
0.00013563176832 


OOTMONMOVAOVVDO DO OCGA DDDD OD OP PPB HHPPYNNNNNWW 


>» 
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» diag (f) 


ans 


—_ 
— 


1.0e+003 * 


» 


4.48533942414282 
4.20598887360221 
3.93601883252608 
2.95655175921535 
2.91232463252813 
2.40665996231507 
2.12840833617064 
2.01547862614741 
1.67501528412810 
1.61012336453201 
1.39849731540630 
1.27385368688073 
1.14741402755860 
1.10615287290734 
0.85542324830144 
0.83293574716167 
0.72369166064969 
0.68994155774074 
0.34608426976888 
0.22304670115122 
0.22787352029346 
0.18559636091107 
0.11394978872364 
0.08854986624844 
0.06796197347123 
0.04772133056046 
0.03761481660823 
0.02913274059718 
0.02217253369177 
0.01562031009177 
0.00771814099641 
0.00652352892696 
0.00285149622307 
0.00178656363092 
0.00046041550163 
0.00016213258818 


Modal Frequencies (with piezos), descending order 


44 














APPENDIX D. MODAL TESTING HARDWARE 


ACCELEROMETER 

1. Model 

Sensotec Model # 3629-05 

2. Available Ranges 

+/- 10g 

3. Sensitivity 

10 mV/g 

4. Usable Frequency (Bandwidth) 
0-250 Hz 

= Excitation 

+/- 5V 

SIGNAL CONDITIONER 

1. Model 

Instruments Division Measurement Group, Inc. Model 2100 
2. Type 

Strain Gage Conditioner and Amplifier System with Integral Power Supply 
3. Gain 

1-2100 

4, Operation 

Full Balanced Operation for Noise Rejection 
IMPULSE FORCE HAMMER 

1. Model 

PCB Piezotronics Model # 086B01 

2. Range 

0-100 Lb. 
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3, Hammer Sensitivity 

52.2 mV/g 

HAMMER POWER SUPPLY 

PCB Piezotronics Model # 483B03 Line Power Supply 

DATA ACQUISITION HARDWARE 

Hewlett-Packard Model HP 3565S Measurement/Data Acquisition System 


46 











LIST OF REFERENCES 


Newman, S., Active Damping Control of a Flexible Space Structure Using 
Piezoelectric Sensors and Actuators, Master's Thesis, Naval Postgraduate School, 
Monterey, California, 1992. 


Feuerstein, M.G., A Comparison of Different Control Methods for Vibration 
Suppression of Flexible Structures Using Piezoelectric Actuators, Master's Thesis, 
Naval Postgraduate School, Monterey, California, 1994. 


Inman, D.J., Vibration: With Control, Measurement, and Stability, Prentice Hall, 
Englewood Cliffs, New Jersey, 1989. 


Meirovitch, L., Elements of Vibration Analysis, McGraw-Hill, Inc., New York, 
New York, 1986. 


Betros, R.S., and Bronowicki, A.J., "Seminar Notes", Active Damping Workshop, 
Spring, 1991. 


Hanagud, S., Obal, M.W., and Calise, A.J.,"Optimal Vibration Control by the Use 
of Piezoceramic Sensors and Actuators", Journal of Guidance, Control, and 
Dynamics, Vol. 15, No. 5, Sept-Oct, 1992, pp. 1199-1206. 


Ewins, D.J., Modal Testing: Theory and Practice, Research Studies Press Ltd., 
John Wiley, New York, New York, 1984. 


Structural Data Research Corporation, /-DEAS™ Test User's Guide, Milford, 
Ohio, 1990. 


James, M.L., Smith, G.M., Wolford, J.C., and Whaley, P.W., Vibration of 


Mechanical and Structural Systems with Microcomputer Applications, Harper and 
Row, Publishers, Inc., New York, New York, 1989. 


~~ 


47 


48 


INITIAL DISTRIBUTION LIST 


Number of Copies 


. Defense Technical Information Center Z 








Cameron Station 
Alexandria, Virginia 22304-6145 


. Library, Code 52 2 
Naval Postgraduate School 
Monterey, California 93943-5104 


. Chairman, Code AA ] 
Department of Aeronautics and Astronautics 

Naval Postgraduate School 

Monterey, California 93943-5002 


. Chairman, Code SP ] 
Department of Aeronautics and Astronautics 

Naval Postgraduate School 

Monterey, California 93943-5002 


. Commander, Naval Space Command l 
Attn: N112 

5280 4th Street 

Dahlgren, Virginia 22448-5300 


. Professor Brij N. Agrawal, Code AA/Ag 2 
Department of Aeronautics and Astronautics 

Naval Postgraduate School 

Monterey, California 93943-5002 


. Professor Hyochoong Bang, Code AA/Ba ] 
Department of Aeronautics and Astronautics 

Naval Postgraduate School 

Monterey, California 93943-5002 


. CDR James C. Trump, USN 1 


6539 Corintia Street. 
Carlsbad, California 92009 


49 


